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Abstract. We consider the evidence for very hard low energy spectra during the prompt phase of Gamma-Ray Bursts (GRB). 
In particular we examine the spectral evolution of GRB 980306 together with the detailed analysis of some other bursts already 
presented in the literature (GRB 911118, GRB 910807, GRB 910927 and GRB 9701 1 1), and check for the significance of their 
hardness (i.e. extremely steep spectral slopes below the EFe peak) by applying different tests. These bursts, detected by the 
Burst And Transient Source Experiment (BATSE) in the ~ 30 keV - 2 MeV energy range, are sufficiently bright to allow 
time resolved spectral studies on time intervals of the order of tenths of a second. We discuss the hard spectra of these bursts 
and their evolution in the context of several non-thermal emission models, which all appear inadequate to account for these 
cases. The extremely hard spectra, which are detected in the early part of the BATSE light curve, are also compared with a 
black body spectral model: the resulting fits are remarkably good, except for an excess at high energies (in several cases) which 
could be simply accounted for by the presence of a supra-thermal component. The findings on the possible thermal character 
of the evolving spectrum and the implications on the GRB physical scenario are considered in the frameworks of photospheric 
models for a fireball which is becoming optically thin, and of Compton drag models, in which the fireball boosts "ambient" 
seed photons by its own bulk motion. Both models, according to simple estimates, appear to be qualitatively and quantitatively 
consistent with the found spectral characteristics, although their possible caveats are discussed. 

Key words, gamma rays: bursts - methods: data analysis - radiation mechanisms: non-thermal, thermal 



> 1 . Introduction 



Since the discovery of these mysterious and puzzling cosmic 
explosions, the time resolved spectral analysis of Gamma-Ray 
Bursts (GRBs) has been a testing ground-floor for the mod- 
els proposed for the 7-ray prompt emission (Ford et al. U995l 
Crider et al. I1997ab . In general GRBs have a continuous cur- 
vature spectrum which is usually characterized (Band et al. 
I1993t by two smoothly joining components defining a peak (in 
a EFe representation of the power per unit logarithmic energy 
interval) at energies around 200 keV. The low energy asymp- 
totic spectrum can be modeled as a power law function in the 
photon distribution N(E) oc E a (i.e. EF E oc E a+2 ), where E 
indicates the photon energy. The high energy spectral compo- 
nent can be described in most cases by a steep power law E° 
(with (3 < a) or an exponential cutoff. The typical values of 
a and (3 are widely distributed around -1 and -2.5 respectively 
(Preece et al. 2000, see also Fig. 1). 

In particular, different authors (Crider et al. 1999 Preece et 
al. 1998a) have stressed the importance of studying the low 
energy component of GRB spectra which is one of the bet- 
ter determined and accessible observables. Crider et al. i 19971 1 
showed that the low energy spectrum evolves in time and that 
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the distribution of the spectral index a extends above - 1 and in 
some cases can be a > 1. 

The most popular radiative process proposed for interpret- 
ing GRB spectra is synchrotron emission by relativistic elec- 
trons in intense magnetic fields (Rees & Meszaros 1994; Katz 
119941 TavaniQ996 ). Optically-thin synchrotron theory predicts 



(Katz 1994) that the low energy photon spectrum can not be 
harder than N(E) oc E~ 2 / 3 . Nonetheless, different authors 
(Crider et al. [T997I Preece et al. I1998al Preece et al. 120021 
have found violations of this limit in many observed spec- 
tra. Following these inconsistencies between the standard syn- 
chrotron model and observations, alternative scenarios have 



been proposed, such as jitter radiation (Medvedev 2001 ); syn- 
chrotron emission from particles with an anisotropic pitch an- 
gle distribution (Lloyd & Petrosian l2000ll2002t : thin/thick syn- 
chrotron emission from a stratified region (Granot, Piran & 
Sari 120001 ): synchrotron self-Compton or inverse Compton off 
photospheric photons (Meszaros & Rees 2000); Compton drag 
(Lazzati et al. 2000 Ghisellini et al. 2000); Comptonization 
of low energy photons by thermal or quasi-thermal particles 
(Liang et al. [T99"7l Ghisellini & Celotti [1999b . The physical pa- 
rameters of these models can be tuned and combined to justify 
the principal observed temporal and spectral characteristics of 
GRBs, producing a considerable number of plausible spectral 
shapes. However, the most extreme (in this case hardest) low 
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energy spectra and their evolution can be used to rule out or 
constrain some of these possibilities. 

Another relevant aspect of GRB spectra is their possible 
thermal character. In fact, the fireball model (Goodman 1986 
Paczynski 1986 1 naturally predicts thermal radiation when the 
fireball becomes transparent. The lack of observational evi- 
dence of thermal spectra motivated the proposal of the inter- 
nal shock model (Rees & Meszaros 1994 > in which the fireball 
energy could be efficiently extracted and transformed into radi- 
ation with a non-thermal spectrum as observed in most bursts. 
Nonetheless, evidence for possible thermal spectra in some 
GRBs has been reported by Preece d200H for GRB 970111 
and Schaefer et al. (poster P-56 presented at the Rome 1998 
Whorkshop on GRBs, private communication). Furthermore 
Blinnikov et al. ([1999 1 propose that also the spectra observed 
in most GRBs could be interpreted as superposition of black 
body spectra with different temperatures. From a theoretical 
point of view thermal emission from the fireball is expected if 
the dominant acceleration agent is internal pressure (Daigne & 
Mochovitch 2002). Alternatively, magnetic acceleration (e.g. 
by Poynting flux) would determine a non-thermal spectrum 
(Drenkhahn & Spruit 120021 but see also Ruffini et al. l2002l . 
The analysis of spectra with a possible thermal character can 
therefore have important implications for the understanding of 
the nature of the fireball acceleration mechanism. 

To the aim of testing the validity/generality of the proposed 
models we have then searched among the 156 GRBs of the pub- 
lished spectral catalog of Preece et al. (2000 1 the bursts with a 
low energy time resolved spectrum harder than N(E) oc E°, 
which represents a limit for most of the above models (see Sec. 
2). We found two extremely hard bursts: GRB 911118 (BATSE 
trigger 1085), whose peak energy evolution has been reported 
in the spectral catalog by Ford et al. ( 1995J, and a new case, 
GRB 980306 (trigger 6630). These two GRBs are carefully 
studied here in terms of their low energy spectral hardness. 
Their spectral evolution on timescales of few tenths of a second 
is presented in Sec. 4.1, 4.2, particularly regarding the slope 
of the low energy spectral component. We also found GRB 
910807 (trigger 647) and GRB 970111 (trigger 5773) whose 
spectral properties have already been presented in the literature 
by Crider et al. Jl997l fT997al l2000l and Frontera et al. {2000) 
(only for GRB 9701 1 1). We includ ed also GRB 910927 (trig- 
ger 829) reported by Crider et al. Jl997l |2000T > which is not 
present in the spectral catalog of Preece et al. (2000 1. These 
bursts have been re-analyzed in order to test the reliability of 
their low energy spectral hardness (Sec. 4.4). In fact, we present 
the tests that we performed in order to determine the statis- 
tical robustness of the spectral hardness, and in particular a 
model independent approach consisting in the comparison of 
each spectrum with a template one of given hardness (Sec. 4.3). 
The evidence that, at least at the beginning of the bursts, the 
spectra are extremely hard suggested their comparison with a 
black body model. This part of the analysis is presented in Sec. 
4.5. The comparison with the low energy spectral limits pre- 
dicted by different models (briefly recalled in Sec. 2 together 
with the previous evidence of hard spectra) is the content of 
the discussion (Sec. 5), where also tentative interpretations of 
the initial quasi-thermal spectral evolution, in the context of 



Table 1. Low energy limiting photon spectral index a [i.e. 
N(E) oc E a ] for various emission models. For clarity the in- 
dices for the energy spectrum F(E) and for its EFe represen- 
tation are also reported. 
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the hot fireball model and of the the Compton drag model, are 
presented. We draw our conclusions in Sec. 6. 

2. Spectral Slopes and Emission Models 

The distribution of the low energy photon spectral index a ob- 
tained from the time resolved spectral analysis of a sample of 
bright BATSE bursts (Ghirlanda, Celotti & Ghisellini l200"2l is 
reported in Fig.[0 The majority of GRB spectra (Preece et al. 
2000 Ghirlanda et al. 2002 ) have a low energy power law spec- 
tral index —3/2 < a < —2/3, i.e. within the limits predicted 
by the optically thin synchrotron model (Fig. |2] dashed lines; 
see also Katz 1994 1, but there is a non negligible fraction of 
bursts (^15%, see also Preece et al. 2000l with a low energy 
spectrum harder than N(E) oc E~ 2 I 3 . The spectral analysis of 
a sample of GRBs by Crider et al. ( 1997 1 also revealed that in 
60% of their cases the spectrum evolves with time. Moreover, 
the spectrum integrated over the pulse rise phase is harder than 
^-2/3 i n 40% of their bursts and can be as hard as TV (i?) oc E 1 
(i.e. EFe oc E 3 ). For comparison the distribution of the a 
parameter for the hard spectra of the bursts considered in this 
work is also reported (gray histogram) in Fig. [2 Note that they 
contribute to extend the distribution toward positive a values 
up to ~ 1.5. 

Let us here just recall that this evidence is hardly rec- 
oncilable with the simplest formulation of the synchrotron 
model and some alternatives have been proposed to account 
for these observations within the frame of this emission pro- 
cess (Papathanassiou 1999 Granot, Piran & Sari 2000 1. Lloyd 
& Petrosian ( 2000 1 propose a scenario in which electrons have 
a small pitch angle distribution (SPD in Fig.|2j, extending the 
range of possible low energy spectral indices produced via syn- 
chrotron to the limit a ~ 0. The same limiting slope can be 
obtained in the "jitter" radiation theory, which is based on syn- 
chrotron emission in a non uniform magnetic field with inho- 
mogeneities on length scales smaller than the electron gyro- 
radius (Medvedev 2001 1. Slopes even harder (i.e. a ~ 1) 
may instead correspond to thermal radiation, such as a portion 
of a black body spectrum (Meszaros & Rees 2000 1 or satu- 
rated Comptonization spectra with Wien peaks which can be 
as hard as a ~ 2 (Liang et al. fl997l Ghisellini & Celotti [T999l 
Ghisellini et al. 120001 . 

A detailed discussion of these models and their comparison 
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Fig. 1. Low energy spectral index distribution from Ghirlanda, 
Celotti & Ghisellini (2002i. The vertical lines represent the 
limits of the emission models reported in Tab. [2 and plotted 
in Fig. |2] The triangles represent the maximum values found 
from the spectral analysis of the bursts presented in this work 
along with their 90% confidence interval (horizontal bars). The 
gray histogram shows the distribution of the low energy spec- 
tral index found from the time resolved spectral analysis of the 
bursts presented in the text. 
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Fig. 2. Low energy spectral limits for different emission models 
(see TabEJ. Examples of the BAND (solid line) and COMP 
model (dot-dashed line) with the low energy photon spectral 
index fixed at -1 [i.e. EFe oc E 1 ] are also represented for 
comparison. 



with our observational findings are presented in sec. 6. As ref- 
erence for the following sections we schematically summarize 
the model predictions in Tab.^ an d display in Fig.|2]their typ- 
ical low energy limiting slopes in a EFe representation. 



3. Data Analysis 

We have analyzed the Large Area Detector (LAD) High Energy 
Resolution Burst (HERB) data which have a high count rate, 
due to the LAD large effective area, and are suited for the spec- 
tral analysis of GRBs for which a moderate energy resolution 
is sufficient to study the typical broad band continuum (Preece 
et al. lT998l . 

We selected the data from the most illuminated detector 
which has the highest signal to noise ratio (S/N), and main- 
tained the instrumental time binning corresponding (in the best 
cases) to single time resolved spectra accumulated for 128 
ms, to follow the spectral variations on the smallest possible 
timescale (thanks to the high S/N there was no need to rebin 
the data before fitting). Each spectrum has been fitted from 
^min ~ 28 keV to -E max ~ 1800 keV which define the typical 
energy window for the LAD data (Preece et al. 1998 2000 1. 

The spectral analysis has been carried out using SOAR 
J1993l > as a quick look tool for the spectral evolution and then 
XSPEC (vll.l) for the individual spectral fitting of time re- 
solved spectra.The background, to be subtracted to each time 
resolved spectrum, has been calculated as the average over a se- 
lected number of background spectra accumulated before and 
after the trigger. 

In the data analysis we have used the standard forward fold- 
ing technique to fit the model spectral functions (convolved 
with the detector response matrix) to the observed spectra. As 
explained in Sec.4, we also used the instrumental spectra in the 
comparison with a simulated spectrum in order to verify the 
robustness of our results. 

3. 1. Fitting models 

We fitted the most commonly adopted GRB spectral models: 
the BAND (Band et al. 1993 1 and COMP model, represented by 
a smoothly connected double power law or by a single power 
law ending with an exponential cutoff, and the sharply con- 
nected double power law model (BPLW). The choice of these 
models is motivated by the fact that they characterize the low 
energy part of the spectrum with a spectral parameter (i.e. the 
spectral index a) which can be simply compared with the pre- 
dictions of the different emission processes. For a detailed de- 
scription of these spectral functions see Preece et al. (2000). 

Due to the extreme hardness of the low energy spectrum of 
these bursts, we decided also to verify if thermal, black body 
like, emission was consistent with the low energy data. We con- 
sidered only the time resolved spectra which have a low energy 
spectral index a > 0.5 when fitted with the BAND or COMP 
models. This indicative value, which is softer than a = 1 (i.e. 
the Rayleigh-Jeans limit of a black body spectrum) has been 
chosen to: a) account for typical errors in the determination of 
a with the BAND or COMP models; b) include the possibility 
that the spectrum could be softer than a = 1, if its peak is at 
low energies, due to the spectral curvature (see also Sec.4. 5). 

The result of each fit was then considered acceptable if the 
reduced x? was lower than 1.5 for typically ^ 110 degrees of 
freedom (dof): in fact the statistical probability of having a bet- 
ter fit is around 0.5 if the reduced Xr is around 1, but a limiting 
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value of Xr ~ 1-5 is suggested if one considers that x 2 depends 
on the quality of the data (Bevington & Robinson 1992). In ad- 
dition we visually inspected the data-to-model ratio to look 
for possible systematics (i.e., sequence of points significantly 
above or below unity). We discarded the black body as a good 
fit when an excess of flux occurred at low energies (i.e. in the 
E 2 part of the spectrum), while we allowed for (a moderate) de- 
viation from unity at high energies (i.e. in the exponential part). 
This latter choice reflects the possible (or even likely) situation 
of having a "supra-thermal" tail to a Maxwellian distribution 
in quasi-thermal plasmas, or a multi-temperature black body 



emission as proposed by Blinnikov et al. ( 1999 



4. Results 

Sec. 4.1 deals with the low energy spectrum of GRB 980306. 
In particular we focus on the early phase of this burst when 
the spectrum is extremely hard and report the time evolution of 
its low energy spectral component. The detailed spectral anal- 
ysis of GRB 911118 (also reported by Ford et al. [T995l onlv in 
terms of the time evolution of the peak energy) is presented in 
Sec. 4.2 with emphasis (again) on the low energy spectral in- 
dex evolution. In Sec. 4.3 we describe the tests we performed 
on the raw data to check for the significance of these results 
on the spectral hardness. In Sec. 4.4 we briefly discuss other 
3 bursts (GRB 910807, GRB 910927 and GRB 970111) with 
hard spectra at low energies which were already presented in 
the literature. We have re-analyzed them, applying the tests de- 
scribed in Sec. 4. 3, to verify their low energy hardness. In par- 
ticular for GRB 910807 we have extended the time analysis 
beyond the first 5 s that were previously studied by Ryde & 
Svensson (1999). The results relative to the black body fits are 
reported separately in Sec. 4.5. 
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4.1. GRB 980306 

GRB 980306 represents a new case of hard burst. It is a single 
peak burst (Fig. [3] top panel) with duration of ~ 6 s, peak flux 
of (17.2 ± 0.4) phot cm~ 2 s _1 at t peak ~ 2.2 s and a total 
(100-300) keV fluence of - 10~ 5 



erg cm 



-2 1 



This burst was 

already present in the Preece et al. (2000) spectral catalogue 
but, as it had been fitted with a BPLW model only, it resulted 
softer than what we have found by fitting its spectra with the 
COMP model. In fact, as already noticed in Ghirlanda et al. 
(2002i, a fit of a model with a sharp break, like the BPLW, to 
a smoothly curved spectrum, although can still give acceptable 
fits, tends to underestimate the low energy spectral hardness. 

We present, for the first time, its hardness evolution, during 
the first 6 s, divided in 14 time resolved spectra. In particular, 
we focus on the extremely hard low energy spectral compo- 
nent. The spectral fits with the different models show that the 
spectra are better represented by the COMP model, since the 
high energy spectrum is very steep (i.e. steeper than a power 
law with spectral index = 10, the limit of the fitting func- 
tion). The time evolution of the low energy spectral index a 
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Fig. 3. Light curve and best fit parameters for GRB 980306. 
Top panel: light curve on the 64 ms timescale, integrated in 
the range 110 - 320 keV and summed over all triggered de- 
tectors. Mid panel: low energy power law spectral index (a). 
Horizontal lines mark the limits a = —2/3 (dot-dashed) and 
a — (dashed). Bottom panel: peak energy -Epoak of the EFe 
spectrum. Error bars represent the 90% confidence intervals on 
the best fit parameters. The filled circles indicate the spectra 
which have been fitted also with a black body model (see Tab. 
Hand Sec. 4.5). 



and of the peak energy E pea k are reported in Fig|3](middle and 
bottom panels, respectively). 

The hardness evolution of a burst is typically described by 
how the spectral index a and i? pca k = {a + 2)Eq change in 
time (here -E pea k is the peak energy of the EFe spectrum and 
Eg is the characteristic energy of the exponential cutoff). The 
evolution of a (mid panel of Fig. [5} indicates that, in the rising 
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phase of the pulse, the maximum hardness (a = 1.1 ± 0.2) is 
reached during the interval t < 1.5 s. The first 4 spectra are 
in fact consistent with a non evolving spectrum (the low en- 
ergy spectral index can be considered equal to its average value 
(a) = 0.73 ± 0.03 at 92% confidence level for the first 1.5 
seconds). After this time the spectrum softens, but it remains 
harder than a = 0. We stress that the low energy spectrum is 
significantly (more than 3a) harder than a = for each of the 
first 7 time resolved spectra (between t = and t = 2.24 s). 
For this specific burst the evolution is hard-to-soft, as shown 
by the peak energy (bottom panel of Fig. [5} decreasing in time 
from ^450 keV to ~50 keV, consistently with the E pea k evo- 
lution reported by Ryde & Svensson ( 2002 1. 

In view of the discussion on the spectral models (sec. 6) we 
also checked the inadequacy of the optically thin synchrotron 
model by fitting the BAND spectrum with the low energy 
power law slope fixed at —2/3. As expected, the residuals at 
low energies show a systematic sign and the reduced Xt is 
larger than 2 indicating that this model is not a good repre- 
sentation of the data. 



4.2. GRB 911118 

The light curve of GRB 911118, reported in Fig.|4](fo/? panel), 
has two main peaks which partially overlap. It is a typical long 
burst with total duration Tgo = 19.2 ± 0.1 s and a background 
subtracted peak flux of (30.6 ± 0.8) phot cm~ 2 s _1 at f pca k = 
6.08 s. Its fluence,in the 100-300 keV energy range, is (27.8 ± 
0.1) x lO^ergcm- 2 . 1 

The spectral analysis that we present covers the first 13 s 
after the trigger time (for a total of 50 spectra), in which the 
time resolved spectrum evolves dramatically with an excursion 



of Aa ~ 1 and AE. D 



400 keV. The best fit is obtained 



with the BAND model. The low energy spectral index is a > 
for the first ~ 6 s (32 spectra, phase A of Fig.0] mid panel) and 
decreases with time (the spectrum softens), but remains harder 
than -2/3 for the time interval up to the end of phase B (see 
Fig. 4). After t ^13 s, the spectrum has a typical a < —2/3 
slope. The hardest spectrum (at t ~ 1.5 s after the trigger) has 
a = 0.74 ± 0.13. The statistical significance of a positive value 
of a for the spectra of phase A is high, being > 3a for most of 
the individual spectra (23/32). The peak energy evolution re- 
ported in Fig.^\(bottom panel) is consistent with what reported 
by Ford et al. (I1995> . According to these fits, the 18 spectra 
belonging to phase B have — 2/3<a<0 and their evolution 
is of the "hard-to-soft" kind, i.e. the peak energy decreases in 
time (becoming smaller than 200 keV in phase B). Fig.|4]shows 
a possible correlation between the peak energy E pea k and a. It 
has been pointed out (Lloyd & Petrosian 2000 1 that such cor- 
relations could result from the effect of the curvature of the 
fitted model when the peak energy is particularly low (e.g. be- 
low ~ 100 keV and close to the low energy spectral threshold). 
However in the case of GRB 91 1 1 18 we can exclude this effect 
because the peak energy is between 500 and 200 keV during 
phase A (Fig|4j. 

Again, we also tried to fit the simplest optically thin syn- 
chrotron model (by fixing the low energy spectral index of the 
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Fig. 4. The light curve and time evolution of the best fit param- 
eters of the BAND model for GRB 911118. Top panel: light 
curve on the 64 ms timescale in the energy range 110 - 320 
keV, summed over all the triggered detectors. The vertical dot- 
ted lines represent the boundaries of the phases in which the 
spectrum is harder than E° (phase A) and E~ 2 / 3 (phase B). 
Mid panel: low energy power law spectral index a, with the 
limits at a = —2/3 (dot-dashed) and a = (dashed). Bottom 
panel: peak energy of the EFe spectrum. The filled circles rep- 
resent the spectra which have been fitted also with a blackbody 
model (see Tab.|2]and Sec. 4.5). 



BAND model to -2/3). In Fig.|5]we report, for comparison, the 
values of the reduced \ 2 r f° r a ^ spectra. The best fit is given 
by the BAND model (asterisks) while the "synchrotron" case 
(triangles) is inadequate to fit these spectra because its \r ^ s 
typically > 2 (for ^110 degrees of freedom) during phase A. 



6 



G. Ghirlanda, A. Celotti & G. Ghisellini: Extremely hard GRBs 



3.0 



2.5 



2.0 



1 .5 



0.5 



A. "A 



A A 
a A A 



A 

/A A 



A 



A 
A 

A A 



4 6 8 

time [sec] 



10 12 



Fig. 5. Reduced x? (f° r typical 110 dof) of the fits with dif- 
ferent models to the spectra of GRB 911118: BAND model 
(asterisks), BAND model with the low energy spectral index 
fixed at -2/3 (triangles). 



The value of Xr f° r the synchrotron model becomes marginally 
acceptable (< 1.5) for most spectra in phase B. 

We conclude that the time resolved spectra of this burst are 
best represented by the BAND model with the 4 fit parameters 
free to vary and that most of the spectra of phase A have a > 
at a > 3fi level. 



4.3. Tests on the Spectral Results 

Considering the relevance of the extreme hardness of these 
bursts for the comparison with the models, we performed dif- 
ferent checks to verify if and at what confidence level they can 
be considered harder than the predictions of the proposed spec- 
tral models (as in Tab. In the following we report on the 
tests performed on the above GRB spectra and their results. 

- reduced fitting energy range: we eliminated from the fit- 
ting energy range the high energy channels to test if they 
substantially influence the low energy fitted spectral slope. 
We restricted the analysis to the range 28-600 ke V (the up- 
per limit being the maximum peak energy found among the 
bursts presented) and compared the results with those ob- 
tained by fitting the spectra on the entire energy window 
(i.e. up to 1800 keV). We found that the best fit spectral 
parameters are consistent within their errors and the values 
of a inferred from the limited energy range case are still 
inconsistent with a = 0. Moreover, Preece et al. ([1998) 
found a soft low energy component in some bursts through 
the analysis of other data type extending below the standard 
LAD threshold of 28 keV. This component could make the 
fitted low energy spectral slope softer than it really is and 
the fit results for the a parameter presented here should be 
considered as lower limits thus reinforcing the evidence of 
extremely hard spectra. We tested this possibility by fitting 
the spectra with a low energy threshold at 50 and 100 keV. 
The comparison with the fits on the entire energy range (28 



- 1800 keV) showed consistent spectral parameters within 
their 90% confidence intervals. 

- background spectrum: if the background spectrum were 
characterized by any feature or absorption edge at low en- 
ergies (e.g. dependent on some electronic failure or anoma- 
lous observational conditions) this might determine partic- 
ularly hard low energy power laws in the background sub- 
tracted spectra. In order to exclude this possibility we com- 
puted several backgrounds selecting different time inter- 
vals, within the same GRB, and found that the background 
spectrum does not depend on the time interval selected for 
its estimate. Moreover, we also fitted the spectra adopting 
different backgrounds calculated for other bursts: the re- 
sults indicate that a is not affected by anomalous features 
in the background spectrum. 

- detector response: the detector response matrix, if incor- 
rectly calculated, could determine a wrong set of best fit 
parameters. We believe this is unlikely as, by fitting the 
spectra with the detector response matrix (DRM) associ- 
ated with another burst but for the same LAD, the results 
are consistent within their 90% confidence level. Moreover, 
although we used the data of the mostly illuminated detec- 
tor, we also tried to fit the data from the other three de- 
tectors that were triggered by these bursts. Although the 
lower signal to noise ratio in these "secondary" detectors 
resulted in a lower statistical significance of the best fit pa- 
rameters, these fits confirmed the extremely hard low en- 
ergy spectral slope of these bursts at least in the second 
best illuminated detector. For the remaining two detectors 
the results are consistent, although the low statistics does 
not allow to confirm the hardness. Note also that Crider et 
al. (1997) showed that combining different detectors can 
make the low energy spectral component softer than it re- 
ally is. 

- comparison with a simulated spectrum: the most direct and 
robust way to verify if the deconvolved spectra are really 
as hard as we found is to compare them with a reference 
spectrum whose shape is well defined. For every time re- 
solved spectrum that we analyzed, we thus simulated a ref- 
erence (template) spectrum with a low energy power law 
slope a ~ and a break energy fixed at the value obtained 
from the fit of the observed spectrum. We divided chan- 
nel by channel the observed spectrum by the template and 
then analyzed the ratio (as long as the observed spectrum 
is harder than the template one their ratio at low energies 
increases with energy). In order to quantify the hardness 
we then fitted such ratio R{E), where E represents the en- 
ergy, with a power-law model R(E) = A E s + B. The 
fitted index 5 is reported in Fig. E0 and Q for GRB 911118 
and GRB 980306, respectively. We note that during phase 
A (the first 6 s after the trigger) of GRB 911118 each of 
the time resolved spectra is at least 3a harder than a flat 
spectrum (i.e. with the same level of confidence that was 
found with the direct fitting method). In the case of GRB 
980306 this test gives a fitted ratio with a poor level of con- 
fidence (< 2a) for each spectrum (with only two spectra at 
> 2a) because of the lower S/N due to count statistics and 
to the propagation of errors in dividing the observed spec- 
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Fig. 6. GRB 91 1 1 18 - Slope of the power law fitted to the ratio 
between the observed time resolved spectra and the template 
spectra (simulated with a = and Ebreak — Eo(fit)). Labels 
A and B refer to the two phases of GRB 911118 (FigBJ and the 
horizontal dotted line is the slope of the reference spectrum. 



trum by the template one. This significance is increased if 
the first 8 spectra are rebinned in time (although the time re- 
solved information is lost) which gives a value of 5 > at 
3.7cr. Furthermore (see also below) also in the case of GRB 
980306 the significance of 6 > (Fig. is higher when 
considering the sequence of S values as a group (since 5 
remains systematically positive). 

As a final consideration let us stress the fact that in the case 
of GRB 911118 (Fig. |4j and GRB 980306 (Fig. EJ there is a 
minority of hard spectra (9 and 4, respectively) which have a > 
at a level of significance lower than 3a. Nonetheless the fact 
that a considerable number of subsequent spectra (32 for GRB 
9 1 1 1 1 8 and 1 2 for GRB 980306) in the first few seconds have 
a spectrum harder (at more than 3a) than E° increases the level 
of significance of the result, if they are considered as a whole. 

4.4. Other hard bursts 

There are other bursts with extremely hard low energy spectra 
for a major part of their evolution, which have already been re- 
ported in the literature. In particular, we found three bursts with 
very hard spectra if fitted by the BAND function. In order to 
have homogeneous results and verify their robustness we have 
then re-analyzed their spectra and applied the tests described 
in Sec. 4. 3. Let us summarize their properties. 

- GRB 910807 has a peak flux of (7.2 ± 0.5) phot cm" 2 
s _1 . The spectrum is harder than a = for the first ~5 s, 
and remains quite constant (~ 0) during the overall burst, 
which lasts about 28 s, as already pointed out by Ryde & 
Svensson ( 1999 1. The peak energy instead shows a tracking 
pattern correlated with the flux and hardens in correspon- 



Fig.7. As Fig.EOfor GRB 980306. 



dence of the second peak (at t — 12 s) as already reported 
by Borgonovo & Ryde (1200 U . 

The complete spectral evolution is reported in Fig. [8] The 
hardest spectrum has a = 1.6 ± 0.5 at 90% confidence 
level. This value is somewhat harder, but consistent with 
what quoted by Crider et al. (2000). A (successful) check 
for the hardness of the first few spectra of this burst has 
been performed using the last method described in the pre- 
vious section. 

- GRB 910927 (Crider et al. l2000lfT997l and GRB 970111 
(Crider et al. 1997a Frontera et al. |2000l are other two 
cases of extremely hard spectra with a low energy power 
law spectral index as high as a = 1.6 ± 0.3. GRB 9701 1 1 
was studied (Frontera et al. 2000) also including the WFC- 
BeppoSAX data. Its extremely hard low energy spectrum 
was confirmed to extend down to 5 keV. Unfortunately, the 
lack of spectral resolution of the GRBM instrument above 
40 keV does not allow to study the complete spectrum, es- 
pecially at its peak. 

In addition we found in the spectral catalogue of Preece et 



al. CPOOl other 3 bursts (BATSE triggers 1974, 2855 and 6350) 
with a > 0, but at a low significance level (< la) and/or for 
a short interval of the light curve. Note that also these cases 
contribute to the a > tail of the spectral index distribution 
reported in Fig.^ 

4.5. Black body spectra 

As anticipated we performed fits on the hardest spectra also 
with a black body function to test if the emission is consistent 
with a thermal model. In fact, the first spectra of all the bursts 
discussed in this work are typically harder than a = 0.5 and 
thus could be consistent with the Rayleigh-Jeans part of a black 
body. Note that the fit could however be unacceptable in terms 
of Xr if the curvature around the peak (i.e. typically in the range 
100 - 300 keV) were broader than a black body one. 
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Fig. 8. GRB 910807 spectral evolution. Same as Fig. [3 



As already mentioned the fits have been performed only 
on spectra with a > 0.5. In order to determine the value of 
this threshold in a we simulated spectra with a black body 
model assuming the detector response of these bursts and a 
typical exposure time ~ 0.2 sec (i.e. the average integration 
time of the spectra presented before). We then fitted these sim- 
ulated thermal spectra with a non-thermal model N(E) oc 
E a exp(— E/Eq) where the power law spectral index a was 
fixed at different values between (a flat photon spectrum) and 
+1 (which best approximates a thermal spectrum). The resid- 
uals of these fits indicated that for a < 0.5 the model spec- 
trum (at low energies) deviates systematically from the data 
(i.e. residuals different from 0) at more than 2a. For this reason 
we selected a = 0.5 as the threshold for fitting the spectra with 
a black body model. 



The criteria applied for the goodness of the fit with the 
black body model were: 

- the reduced \r should be smaller than a fixed reference 
value of 1.5. 

- the data-to-model ratio can be systematically and signifi- 
cantly greater than 1 only at energies above the peak, as a 
possible consequence of the presence of a supra-thermal or 
multi-temperature spectral component (whose nature is not 
discussed here). 

The results of the fits for each spectrum are reported in Tab. 
13 the start and stop time, with respect to the trigger time, are 
given in col. 2 and 3. The model temperature kT (keV) and its 
normalization (erg cm -2 s" 1 ) are reported in col. 4 and 5. The 
value of the reduced Xr /dof with the corresponding degrees of 
freedom and the photon flux computed from the fitted model in 
the 28-1 800 keV energy range are given in col. 6 and 7, respec- 
tively. Column 8 and 9 report the spectral index (a) obtained by 
the fit with the non thermal (BAND or COMP) models and the 
associated Xr/dof. 

For illustration, we report in Fig.|9]some examples of fitted 
spectra for GRB 980306 and the corresponding data-to-model 
ratios. The first two spectra are considered well fitted by the 
model because Xr = 1-19 and 0.9 for 111 dof, respectively, 
and the ratio is systematically greater than 1 only at high ener- 
gies. In the other two cases the black body fits give less accept- 
able results because the data-to-model ratio shows systematic 
deviations from 1 at low energies. 

In most of the bursts the best fits with the black body model 
are obtained at the beginning: in the later stages (e.g. after the 
first 2.5 s in the case of GRB 980306) the low energy com- 
ponent softens and becomes incompatible with the black body 
model as the reduced Xr an d the data-to-model ratio indicate. 

The black body temperature and total flux (as derived by 
the model) evolve in time as reported in Fig. ^3 f° r a H the 
good black body spectral fits. Note the decrease of the temper- 
ature with time in almost all these bursts, which in the figure 
are compared with a dependence Tbb °c i~ 1//4 (dotted line). 
Although the temperature decreases, the flux slightly increases 
or remains constant in the first phases, and rapidly decreases 
thereafter. 

Before concluding we note that although not publicly avail- 
able, other possible evidence of a black body emission from 
GRBs has been reported in poster proceedings (Palmer D., pri- 
vate communication). The possible thermal black body spec- 
trum of GRB9701 1 1 has also been recently presented by Preece 
(2001 ): the spectrum integrated over its first 5 s could be fitted 
by a black body with kT = 55 keV which is consistent with 
the average value of the black body temperatures that we ob- 
tain from the fits of the time resolved spectra covering the same 
time interval CFigllOi. 

5. Discussion 

GRB 980306 and 91 11 18, together with GRB9 10807, 910927 
and 970111 already reported in the literature, bring new con- 
straints on the extremely hard spectra that the emission process 
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Fig. 9. Examples of black body fits for GRB 980306 (see Tab|2j. The integration times are indicated for each plot. The data-to- 
model ratios are reported in the bottom panels. 



responsible for the prompt phase must be able to satisfy. These 
bursts represent the cases with the hardest low energy power 
law spectral index and significantly contribute in extending the 
a distribution reported by Preece at al. ( 2000 1 and Ghirlanda et 
al. d2002> to values greater than 0, up to a ~ 1.5. For the bursts 
studied here a, i.e. the spectral index of the low energy photon 
spectrum, remains positive for a major part of the pulse/s, like 
in the case of GRB 91 1 1 18 or GRB 9701 1 1, or for its first ris- 
ing part (GRB 980306 and GRB 910807). 

In the following we consider several emission models pro- 
posed so far to see if they can account for these hard spectra, 
except for the photospheric and Compton drag model which 



will be discussed in Sec. 5.5.2 in relation to the quasi-thermal 
character of the initial phases of these extremely hard bursts. 



5. 1. Synchrotron emission 

Optically thin synchrotron emission — The optically thin 
synchrotron model, in its simplest formulation, can produce 
at most a spectrum as hard as a = —2/3 (Katz 1994 1. 
Furthermore if one considers the very short cooling timescales 
(much shorter than any conceivable current exposure time) 
the predicted 'cooled' spectrum has a = —3/2 (see also 
Ghisellini, Celotti & Lazzati 2000). 
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Fig. 10. Total flux (top panel) and black body temperature (bot- 
tom panel) as a function of time for the spectra reported in Tab. 
[2]and considered acceptable in terms of residuals. The dashed 
lines in the bottom panel correspond to Tbb oc t^ 1 / 4 . 



Synchrotron self-absorption — One alternative, investigated 
among others by Papathanassiou ( 1999 ), is that the medium re- 
sponsible for the synchrotron emission also absorbs these pho- 
tons. Also Granot Piran & Sari (2000) have proposed a model 
based on synchrotron self absorption by a stratified electron 
population. In these cases it is possible to produce low energy 
spectral slopes as hard as a = 1 (for a thermal or a power 
law electron energy distribution with a low energy cut-off) or 
a = 1.5 (for a power law extending to low energies). 

However, for synchrotron self-absorption to be effective 
in the typical BATSE energy range, a very high density of 
relativistic electrons is necessary. Assuming that the electron 
distribution extends to low energies (this minimizes the re- 
quired particle number) with a power law distribution of in- 
dex p (i.e. N c (j) oc 7~ p , where 7 is the Lorentz factor of 
the electrons), the synchrotron self-absorption frequency v t in 
the comoving frame, assuming p = 2, is (e.g. Ghisellini & 
Svensson 1991) v t ~ 2 x 10 5 ^ B {tt/B) 1 / 3 Hz, where t t 
is the Thomson optical depth and v-q = eB / (2irm e c) is the 
Larmor frequency. In order to have v t ~ 10 16 Hz, correspond- 
ing to ~ Tv t in the BATSE range, t t B 2 ~ 5.7 x 10 12 ^ t 3 16 . 
The Comptonization parameter y ~ tt(7 2 } corresponding to 
the required particle density, would largely exceed unity (mak- 
ing inverse Compton process more efficient than synchrotron) 
unless the magnetic field is extremely high (B ~ 10 7 G). 
This would in turn imply a (magnetically dominated) fireball 



with isotropic (equivalent) energy in excess of 10 54 erg, us- 
ing E ~ ATrR 2 (B 2 /8TT)r 2 ct burst and assuming T ~ 10 2 and 
B ~ 10 13 cm. A magnetic field of ~ 10 7 G has been derived 



also by Crider & Liang ( 1999) from the fit of a simplified self 
absorbed synchrotron model to the spectrum of GRB 9701 11. 



Small pitch angles — Lloyd & Petrosian (2000) have pro- 
posed that in rarefied, highly magnetized plasmas, turbulence 
can give origin to an anisotropic pitch angle distribution of 
emitting particles. In this case the optically thin synchrotron 
spectrum is modified at low energies (compared to the stan- 
dard F(E) oc £ 1/3 ), with a limiting slope F(E) oc E. This 
can accommodate some, but not all, of the hard spectra that 
we presented. Note also that according to this scenario to pre- 
vent strong cooling the particles have to be re-accelerated, con- 
trary to one of the basic assumptions of the internal shock sce- 
nario, and in order not to be re-isotropized by scattering, the 
inverse Compton process must be much less efficient than the 
synchrotron one. 

5.2. Jitter radiation 

A variant of the standard synchrotron emission theory, the jitter 
radiation (Medvedev 2001 1, can justify a flat low energy spec- 
tral slope. It originates as the emission of relativistic electrons 
in a non uniform magnetic field with inhomogeneity length 
scale smaller than the Larmor gyro-radius. Medvedev (2001) 
proposed a composite model for GRB spectra by assuming the 
presence of a small scale magnetic field, which causes the jit- 
ter radiation component, and a larger scale field producing a 
standard synchrotron spectrum. The composite spectrum has a 
broad bump in correspondence of the jitter radiation character- 
istic energy, which depends only on the magnetic field proper- 
ties and not on the electron energy and if the small scale field 
prevails, the spectrum has a limiting a = 0. This implies that 
only a minority of the hard spectra reported here could be con- 
sistent with such a model. 



5.3. Compton attenuation 



Compton attenuation (Brainerd et al. 1998) can produce a low 
energy (20-100 keV) hard component because of the propaga- 
tion and scattering of the intrinsic spectrum by material (not 
participating to the bulk flow) located between the fireball and 
the observer. The original high energy spectral shape can be un- 
altered because of the decline of the Klein-Nishina cross sec- 
tion with energy, while at lower energies, where the scattering 
process is more efficient, the spectrum is hardened. The most 
appealing features of this model are that it could explain the 
clustering of the peak energy of bursts around some hundred 
keV and produce a flat slope at low energies. But in order for 
this to work, the Thomson optical depth of the scattering ma- 
terial must be quite large (of the order of tt ~ 10) and this 
would smooth out the observed burst light curve by smearing 
the variability on the smallest time scales. Furthermore, by de- 
creasing the transmitted (spectrally unaltered) flux, this mech- 
anism requires intrinsic powers greater than what we observe, 
exacerbating the energy budget problem. 

5.4. Quasi-thermal Comptonization 

The above models are based on the assumption that the emitting 
particles are highly relativistic, as a result of "instantaneous" 
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acceleration at a shock front. However, if the heating mecha- 
nism operates on timescales of the order of the dynamical time 
(or the light crossing time of a shell in the internal shock sce- 
nario), then the heating and cooling processes could balance, 
leading to typical energies of the electrons that are sub or only 
mildly relativistic (Ghisellini & Celotti 1999). In this case the 
dominant radiation process is Comptonization of seed soft pho- 
tons by a quasi-thermal particle distribution. 

Different models assume different sources for the seed pho- 
tons and different parameters (i.e. typical size for the emitting 
region, magnetic field, role of the photospheric radiation, and 
so on) within the internal shock scenario. 

Liang et al. (1997), among the first to propose such process 
for the prompt emission, considered a relatively large emit- 
ting region embedded in a relatively weak magnetic field. This 
model was applied to GRB 990123 by Liang (1999), who found 
quite extreme best fit values for the required total power and 
number of emitting particles. Ghisellini & Celotti (1999) pro- 
posed that Comptonization acts on the self-absorbed cyclo- 
synchrotron radiation produced by the same quasi-thermal 
particles responsible for the multiple scattering process, and 
pointed out the role that electron-positron pairs can have in 
keeping the temperature (or the mean energy) of the particles 
within a narrow range. On the same line Meszaros & Rees 
(2000) pointed out the importance of the residual photospheric 
emission as source of soft photons. 

The basic features of the Comptonization model is that, 
in the quasi-saturated regime a Wien peak with its character- 
istic F(E) oc E 3 shape can form at high energies (i.e. for 
hv ~ kT). In principle this can thus explain very hard spec- 
tra. One possible difficulty of this model is the slope at lower 
energies, whose saturated value should be a = —1. A further 
possible problem for quasi-thermal Componization is that the 
temperature of the emitting particles (in the comoving frame) 
is expected to be of order of ~ 50 keV, which would lead to 
an observed E peak of few MeV (Ghisellini & Celotti 1999). 
This is a somewhat high value, also in the light of the results 
presented in this work, which show that even spectra time re- 
solved on scales < 1 sec present a peak energy of E pea k < 500 
keV. 



5.5. The possible thermal character of the initial 
phase 

The first 1-5 s of the emission of the bursts presented in this 
work were found to be consistent with a black body spectrum 
with typical temperatures initially around ~ 100 keV and de- 
creasing to 30 - 40 keV. After this initial phase a softer (non- 
thermal or multi-temperature) character of the spectrum be- 
comes dominant. Clearly this behavior can be a powerful di- 
agnostic. In particular it seems to favour two among the pro- 
posed models, namely photospheric and Compton drag emis- 
sion, which both predict - as recalled in Sec. 5.5, 5.6 - a ther- 
mal spectrum during the first phases. Here we briefly discuss 
which physical constraints can be quantitatively gathered from 
the data. 



5.5.1. Photospheric emission 

When the fireball is becoming optically thin, during the accel- 
eration or coasting phase, its internal energy is emitted with 
a black body spectrum whose observed temperature is blue- 
shifted by a factor T with respect to the comoving one. The pos- 
sible presence of a photospheric component in the spectrum of 
GRBs has been theoretically investigated by Meszaros & Rees 
(2000 1. Recently, Daigne & Mochkovitch (2002) have recon- 
sidered this possibility, finding quite tight limits on the model 
of hot fireballs accelerated through internal radiation pressure, 
imposed by the absence of an initial emission phase (or precur- 
sor) with a black body shape. 

In this scenario we have two options: we can relate the ini- 
tial phase to a single fireball, in the process of becoming opti- 
cally thin, or - following Daigne & Mochkovitch (2002) - to 
an ensemble of Af shells, each becoming optically thin at ap- 
proximately the same distance. In both cases at transparency: 

Ef or 



2n0 2 Rh 



2 TN 



1 



(1) 



where R t is the transparency radius of the fireball, collimated 
into two cones of semi-aperture angle 0, E{ is the total fireball 
energy (a fraction e 7 of which is radiated as photons), and <tt 
is the Thomson cross section. Assuming that the surface of the 
fireball is emitting, at R t , as a black body with a comoving 
temperature T' = T b s /r, we have 
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where (.Lbb) is the average black body luminosity in the ob- 
server frame for the time £bb and a is the Stefan-Boltzmann 
constant. Both eq. (2) and (3) show a Ef oc 8 2 R 2 dependence, 
allowing to solve for the bulk Lorentz factor: 
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with X^s = 10 9 T o bs,9 K (corresponding to the observed black 
body peak of ~ 100 keV). We stress that this estimate of T is 
independent of the degree of collimation of the fireball. 

In the case of a single shell t^B is the time needed for the 
fireball to become transparent, implying 
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Therefore, if the emission comes from a single shell, the trans- 
parency radius is very large, implying unreasonably large val- 
ues for the luminosity and energy, i.e. from eq.(3), Ef ~ 
3 x 10 58 £gg 3 &7 5 ^ 3 erg, where 9 is expressed in 

units of 0.1 radians. 

If, instead, the observed emission comes from a series of 
shells the transparency radius can be much smaller, allow- 
ing a "reasonable" luminosity of the photosphere: for Ef — 
10 51 i?f ,51 erg we have 
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We conclude that a series of shells (of total energy Et) each 
becoming transparent at R t can account for the observed black 
body emission. The case of J\f shells for large JV is clearly 
equivalent to a continuous or quasi-continuous flow. 

As Daigne & Mochkovitch (2002) pointed out, and in 
agreement with our findings, the photospheric radiation is 
likely to be visible only during the first phases of the burst light 
curve, as long as the optical depth of the material ahead of the 
shell which is releasing its thermal radiation is negligible and 
before internal shocks take over. 

5.5.2. Compton drag 

If the circum-burst environment is characterized by quite a 
large photon density, as is the case of bursts following a su- 
pernova explosion, or for fireballs produced in the matter- 
evacuated funnel of an hypernova, there can be a strong in- 
teraction between these seed photons and the fireball itself, as 
postulated in the so-called Compton drag model. The ambient 
photon energy is boosted by the factor T 2 at the expense of the 
fireball kinetic energy (Lazzati et al. 2000 1. If the funnel or the 
young supernova remnant are characterized by a single temper- 
ature Tsn, and if the fireball does not decelerate, the emitted 
spectrum is a black body at a temperature ~ 2T 2 Tsn- If in- 
stead the seed photons have a range of temperatures (as likely 
to be in the case of a funnel, hotter in the central parts), and/or 
the process is so efficient to decelerate the fireball, than the fi- 
nal spectrum will be a superposition of the locally produced 
black body spectra, as calculated by Ghisellini et al. (2000 1. 
The resulting spectrum can indeed resemble the very hard spec- 
tra found, at least for the initial part of these bursts. In fact, once 
the first shells have swept up the seed photons, the efficiency 
of Compton drag is greatly reduced, since the time to replen- 
ish the circum-burst environment with seed photons is typi- 
cally longer than the duration of the burst itself. On the other 
hand the Compton drag model favours the formation of internal 
shocks, since after the first shells deceleration, the subsequent 
ejecta can more easily collide with them and produce shocks. 
Then another radiation source can become efficient after the 
first phase (lasting ~1 to a few seconds). This possibly also 
explains the hard to soft evolution: at first the spectrum, due 
to Compton drag, is very hard at low energies (F(E) oc E 2 ), 
whereas at later times the spectrum produced in internal shocks 
through other processes becomes softer (F(E) oc E 1 ^ 3 or 
oc E°). 

Let us then assume that the time interval for which the black 
body lasts corresponds to the emission from a single shell and 
that the emission peaks when the shell becomes transparent (eq. 
2), implying R t ~ ct t T 2 . 

The total energy Ecn due to the Compton drag process, 
produced in the observed time t t , is of the order 

E CB ~ 2n9 2 RU2T 2 )a(^y 

~ 3 x 10 51 ^ lt 3T o 4 bS)9 erg (6) 

where a — 7.56 x 10 _15 erg cm -3 K~ 4 is the radiation constant. 
In this form Ecn does not depend on T (although of course 



T$n does). We can then derive the transparency radius (from 
eq.(2) and (8)) and the bulk Lorentz factor: 

l 8e 7 TO p J 
~ 3xl0 14 r o 8 £ 9 4 /5 S 2/5 cm (7) 

r= (f) 1/2 -i00T o t 9 4 /5 s 1/5 (8) 

for a temperature of the seed photons corresponding to Tg^ ~ 
5 x 10 4 K. 

These simple (and rough) estimates indicate that the 
Compton drag model is viable. 

5.5.3. Time evolution 

Although the limited number of time resolved spectra and 
bursts do not allow to search for a general evolutionary be- 
haviour, we can briefly comment on the spectral and dynamical 
evolution of the initial thermal phase of these 5 bursts. As re- 
ported in Fig.^Sin a ll the 5 bursts, the black body temperature 
decreases with time (bottom panel) and its evolution is intrigu- 
ingly close to T b s oc t~ x l A for the first few seconds, and then 
drops. In the same time interval the luminosity (Lbb) remains 
constant or mildly increases (top panel of Fig. llOt . We therefore 
observe, at least in the first phase, a decrease in the observed 
temperature without a corresponding decrease in the observed 
black body flux. We plan to examine more deeply these obser- 
vational behaviour in a future work: here we just briefly com- 
ment on it in relation to the scenarios discussed above. 

In the context of the photospheric model discussed in Sec. 
5.7.1, the black body evolution could be explained by succes- 
sive shells having an increasing baryon loading and smaller 
Lorentz factor, thus becoming transparent at increasing dis- 
tances R t . This can cause the observed temperature to decrease 
(because the Lorentz factor is smaller) without a decrease in 
flux (because the radius is larger). 

In the context of the Compton drag model, a decrease in the 
observed temperature can be due to the deceleration of the fire- 
ball and/or a decreasing temperature of the seed photons with 
distance. If we consider the luminosity due to the Compton 
drag process (Eq. (8)), we have Lcd °c Eco/t oc t 2 T^ hB , 
which predicts Lqd oc £ if the time behavior of the observed 
temperature is indeed T D b s oc i -1 / 4 . A weaker dependence of 
the observed flux on time can occur if the fireball is becoming 
transparent (i.e. only a fraction tt of the seed photons can be 
scattered) with tt decreasing with distance and/or if the fireball 
distance becomes larger than the typical dimension occupied 
by the seed photons. 

6. Conclusions 

We have presented the spectral evolution of GRB 911118 and 
of a new case of hard burst, GRB 980306, together with some 
other hard bursts already reported in the literature. Their low 
energy spectral component is harder than N(E) oc E° for a 
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considerable period of their main peak emission and even after 
this initial phase their spectrum remains harder than the syn- 
chrotron limit E~ 2 I Z . We have applied different tests to verify 
the significance of the low energy hardness, concluding that 
these results are indeed robust. 

These GRB prompt spectra represent a challenge for the 
proposed emission scenarios as shown through a comparison 
with the limiting spectral shapes predicted by such models. We 
pointed out the difficulties that synchrotron emission, even in- 
cluding the effects of self-absorption, small pitch angles parti- 
cle distributions and jittering, has in explaining spectra harder 
than N(E) oc E°. Comptonization models also have difficul- 
ties, even if they are consistent with very hard spectra in a lim- 
ited range of energies. Of course, such a pruning of the forest 
of emission scenarios refers only to the phases characterized by 
the harder emission, while it is possible and likely that different 
processes dominate at different stages of the prompt evolution 
and/or in different GRBs (see also below). 

The main result of this study is the possible thermal char- 
acter of the first emission phase of all the bursts we consid- 
ered and the findings on the evolution of such a thermal phase. 
The conventional scenario, of a fireball accelerated by its own 
internal pressure, indeed predicts such an initial thermal char- 
acter when the fireball becomes transparent. It was indeed its 
absence in previously considered bursts that led Daigne & 
Mochkovitch (2002) to favour a rather cold fireball scenario, 
where at least part of the acceleration was due to magnetic 
forces. In our bursts, on the contrary, the luminosity in the 
thermal phase is a significant fraction of the total, therefore 
in agreement with the hot fireball scenario. However, the fire- 
ball could be cold initially, and be heated later (but before it 
becomes transparent). The heating agent could be magnetic 
reconnection (see e.g. Drenkhahn & Spruit 2002) or internal 
shocks occurring at early phases. The latter case demands that 
the typical bulk Lorentz factors of the shells are small enough 
for the shells to collide before the transparency radius. 

Alternatively, a cold fireball could work if the thermal spec- 
trum we see is produced by Compton drag, in which circum- 
burst radiation is boosted to high energies by the fireball bulk 
motion. 

At present, we are not able to discriminate between the pho- 
tospheric and the Compton drag scenario. However, a key dif- 
ference between the two scenarios is the fact that the seeds pho- 
tons for the Compton drag process can be "used" only by the 
first shells, because the time needed to refill the scattering zone 
with new seeds exceeds the duration of the burst. Therefore 
observing blackbody emission for a long time, or during the 
rising phase of two time resolved peaks would be difficult to 
explain in terms of the Compton drag process. We finally note 
that emission at times greater than a few seconds can well be 
due to other processes, possibly linked to internal shocks start- 
ing to dominate at later phases. 
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Table 2. Black body fits. 
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